Implications of ultra-high energy neutrino flux constraints 
for Lorentz-invariance violating cosmogenic neutrinos 
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We consider the implications of Lorentz-invariance violation (LIV) on cosmogenic neutrino observations, 
with particular focus on the constraints imposed on several well-developed models for ultra-high energy cos- 
mogenic neutrino production by recent results from the Antarctic Impulsive Transient Antenna (ANITA) long- 
duration balloon payload, and Radio Ice Cherenkov Experiment (RICE) at the South Pole. Under a scenario 
proposed originally by Coleman and Glashow, each lepton family may attain maximum velocities that can ex- 
ceed c, leading to energy-loss through several interaction channels during propagation. We show that future 
observations of cosmogenic neutrinos will provide by far the most stringent limit on LIV in the neutrino sector. 
We derive the implied level of LIV required to suppress observation of predicted fluxes from several mainstream 
cosmogenic neutrino models, and specifically those recently constrained by the ANITA and RICE experiments. 
We simulate via detailed Monte Carlo code the propagation of cosmogenic neutrino fluxes in the presence of 
LIV-induced energy losses. We show that this process produces several detectable effects in the resulting atten- 
uated neutrino spectra, even at LIV-induced neutrino superluminality of (« v — c)/c ~ 10~ 26 , about 13 orders of 
magnitude below current bounds. 



In the current Standard Model (SM) of particle physics, 
neutrinos are massless and unmixed, and couple to other fun- 
damental particles only through the weak interaction. This 
theory is very successful at describing the gauge boson fields 
and their interactions, and the three generations of quarks and 
leptons. However, observations of neutrino oscillations imply 
mixed neutrinos with non-zero mass, the most compelling ev- 
idence requiring extensions to the SM. The search for a funda- 
mental theory that can explain and predict properties of neutri- 
nos with non-zero mass has in turn become one of the most ac- 
tive areas of high energy particle physics. It is thus natural to 
anticipate that neutrino observations may provide unique op- 
portunities for investigation of other beyond-standard-model 
physics such as Lorentz invariance violation (LIV). While a 
recent claim of LIV in the form of superluminal behavior in 
muon neutrinos [ 1 1 appears to have been spurious and due to 
subtle instrumental effects, the possibility that LIV can ap- 
pear at very small levels in many sectors of particle physics 
has been seriously explored for decades. Investigation of LIV 
that could lead to non-tachyonic superluminal particle veloc- 
ities has received renewed motivation from recent efforts to 
develop consistent quantum gravity theories, through explo- 
ration of the possibility that Lorentz invariance may not be an 
exact vacuum symmetry at high energies EJEl. 

SM modifications which specifically address LIV have 
been developed both as a way of incorporating the low-energy 



effects of spontaneous CPT-violation into the SM [4], and for 
high-energies, via general perturbation analysis [5, 6|. The 
former work has been generalized to a Standard Model Exten- 
sion (SME), and as such it provides a comprehensive frame- 
work for analysis of almost any LIV effect in any particle 
sector; recent compilations of experimental constraints now 
give limits for hundreds of different possible LIV parame- 
ters ||8][10), and particular applications to the neutrino sector, 
with a primary focus on low-energy phenomena, have also 
been developed within this framework Q, and more recently 
for even a fully relativistic formulation with operators of arbi- 
trary dimension ||9] . 

The work of Coleman and Glashow |UJ |6), which fo- 
cused on the high-energy limit, provides some specific predic- 
tions regarding energy-loss mechanisms for superluminal LIV 
charged particles. These results have been used to provide 
stringent constraints on LIV on both ultra-high energy cosmic 
ray (UHECR) protons and UHECR photons [11]. More 
recently, similar considerations have been extended to cosmo- 
genic neutrinos arising in intergalactic UHECR propagation 
showing that LIV in the hadron sector could lead to suppres- 
sion of photomeson production, and a resulting decrease in 
the daughter neutrino fluxes lfT2ll . Thus detailed UHE cosmo- 
genic neutrino measurements and constraints can indirectly 
signal very small levels of LIV in the hadronic sector. 

In the neutrino sector, LIV can appear in a very large va- 
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riety of forms; for example, the recent comprehensive study 
by Kostelecky & Mewes gives 369 possible coefficients of 
LIV parameters for models which are restricted to be renor- 
malizable (9). In general, these coefficients can be associated 
with both flavor-changing effects, where the velocity between 
neutrino mass eigenstates can differ, and flavor-blind effects, 
where all flavors show an LIV-induced effect. In addition, di- 
rectional dependence of the particle velocity may be implied 
when the LIV involves a preferred inertial frame. The flavor- 
dependent effects may appear as a type of oscillation signal 
for neutrinos propagating along different pathlengths or along 
different directions, but need not imply superluminal LIV. De- 
tection of an LIV-induced signal due to these effects requires 
a search for modulation of the neutrino signal as a function of 
of the effective propagation distance, and such searches can 
be extremely sensitive since they search for differential ef- 
fects. In contrast, flavor-blind effects would be responsible 
for producing LIV that would be unique to the entire neutrino 
sector, affecting each mass eigenstate in the same way, and 
potentially leading to superluminal states. Searches for these 
effects must differentiate the neutrino sector velocity against 
another sector such as the photon sector, and are thus much 
more difficult in practice. 

Currently, velocity differences between different neutrino 
mass eigenstates are constrained to parts in 10~ 23 by MI- 
NOS ifTSl with accelerator-based neutrinos, and to parts in 
10~ 27 by IceCube using atmospheric muon neutrinos. How- 
ever, constraints on the overall velocity difference of all neu- 
trino flavors with respect to photons are many orders of 
magnitude weaker, both than the constraints on intra-mass- 
eigenstate variation, and compared to constraints on many 
other particle sectors in the SME. 

In this work we focus only on LIV which can lead to an 
overall LIV superluminal behavior for the neutrino sector, and 
we will concentrate only on isotropic effects without direc- 
tion dependence. The Lagrangian density prescribed in the 
LIV-motivated SME leads to a modified energy-momentum 
relation for neutrinos J6] Q3] : 

E 2 = m 2 v c 4 + [l+28 v ( Pv )}c 2 p 2 , (1) 

where E v ,m v ,p v are the particle's energy, mass, and momen- 
tum and the parameter Sy(pv), which is in general a spin- 
dependent, linear combination of several SME parameters, 
and can be identified with the dimension-three isotropic co- 
efficient in the Standard-Model Extension denoted (cl)oq, or 
(c /)oo in the flavor-blind, oscillation-free limit ll9l [T6l . 

This term effectively determines the maximum possible ve- 
locity in direction /3 V as u v w c(l + §v), where we have kept 
only terms to first order in 8 V . This relation implies kinemati- 
cally that the velocity of a neutrino may exceed that of photons 
in vacuo in the high-energy limit. This behavior is distinct 
from the hypothesis of Lorentz-invariant tachyonic neutrinos, 
first proposed by Chodos et al. ifTTll . which have very different 
phenomenology and constraints; see 1 18 1 for a recent review. 

Recent work |[T9l has shown that LIV-induced superluminal 
neutrino propagation in vacuum will lead to several energy- 
loss mechanisms not present in standard model propagation. 
One mechanism in particular, electron-positron pair creation, 



or pair bremsstrahlung, is efficient enough for even very small 
values of 8 V , that it leads to strong constraints in any exper- 
iment where neutrinos of sufficient energy are detected af- 
ter propagation over significant distances. This in turn im- 
plies that even miniscule non-zero values of 8 V will effectively 
attenuate cosmogenic neutrino fluxes to undetectable levels. 
Detection of such effects take the form of a disappearance 
search, coupled with the possibility of modification of the par- 
ent neutrino spectral energy distribution. This then alleviates 
the need to compare the neutrinos against another particle sec- 
tor to establish the propagation characteristics. 

Cowsik et al. [20 1 have recently investigated the bounds 
on LIV in the neutrino sector that arise from atmospheric 
muon neutrinos observed in IceCube, using the convention 
E v = p v c(l + a v ) for the LIV parameter, where they have ne- 
glected the neutrino mass in the high energy limit, and thus 
a v = 8v as used above in the SME. At high energies in the 
pair bremsstrahlung process noted above, V; — > vj + e~ + e + , 
the neutrino loses about 3/4 of its energy, and thus to first order 
a single pair bremsstrahlung interaction during neutrino prop- 
agation may be treated as an effective neutrino decay |[T9ll20ll 
with characteristic decay time 

T v = TcG £ v ~GeV ^ S - ( 2 ) 

where %G = 6.5 x 10~ n s. For IceCube's observations of 
upcoming atmospheric neutrinos, by requiring that the decay 
time exceed the crossing time of the Earth, Cowsik et al. find 
a resulting bound of oc v < 1CT 13 . 

Cowsik et al. also use similar analysis to derive a first- 
order estimate of the possible upper limit to a v that would 
arise from any observation of UHE cosmogenic neutrinos, by 
again requiring that the effective decay time be longer than 
the propagation time from the cosmogenic neutrino produc- 
tion site to Earth. However, they considered only neutrinos 
arising from sources within the nearest 100 Mpc. While this 
is the most conservative approach, it is in tension with the 
typical scenario presented by most UHE cosmogenic neutrino 
models. In fact, in all current models for their production and 
propagation, only a negligible fraction of cosmogenic neutri- 
nos observed at Earth arise from nearby sources; by far the 
largest contribution is from relatively high-redshift sources. 
To make an estimate consistent with this fact, we must first 
determine how to treat the neutrino propagation and possible 
effective decay over cosmological distances. 

At a given redshift z, the propagation time for photons to 
the current epoch is given by the lookback time tf, 

k{z) = f\(z')dz' (3) 
Jo 

where r|(z) =dt/dz = t H [(l + z)'£(z')]~ 1 |2T). Here t H = H Q 1 
is the Hubble time (t H ~ 3.09 x I0 ll h sec for H = 100ft km 
s Mpc -1 , h ~ 0.7) and the function £(z) is given by 11221 : 

<E(z) = [ClM(l+z) 3 +Q. k (l+z) 2 +Q.A} 1/2 . (4) 

Here the curvature parameter £2^ = 1 — Q.m — £2a, where 
Q.m,Qa are the matter density and cosmological constant, re- 
spectively. We consider three possible cosmologies, Einstein- 
deSitter (£2m = L^a = 0), a low-density, high-curvature 



3 




high-lambda 

low-density 

Einstein-deSit 



0.05 0.1 



0.5 1 
redshift, 



10 



FIG. 1: Propagation time as a function of redshift z for neutrinos 
for 3 different cosmological models. The dotted line also indicates 
Euclidean propagation time. 



model (£2^ = 0.05, £2a = 0), and the currently favored high- 
lambda model (Q.m = 0.3, Ha = 0.7). Fig. [T] shows the look- 
back time for these three models, for z > 0.03, approximately 
130 Mpc or more. The dotted line also shows the Euclidean- 
space approximation, t =ztn- At z > 0.1 departure from 
Euclidean propagation is evident, and at z ~ 1 the propaga- 
tion time begins to saturate asymptotically to the Hubble time, 
with less than a factor of two variation among the different 
cosmological models. 

A large number of studies of UHE cosmogenic neutrinos 
have been published since 1969 when Berezinsky and Zat- 
sepin [23] first described their origin as secondaries in the 
Greisen-Zatsepin-Kuzmin (GZK 11241 [25 1) process by which 
UHE cosmic ray protons resonantly scatter off the cosmic 
microwave background photons. For our purposes here, we 
refer only to several recent results that provide estimates 
of cosmogenic neutrino spectra as observed at Earth, for a 
range of model parameters, and for sources of different red- 
shift [EHUD. 

Kotera et al. 11261 present a decomposition of the neutrino 
energy spectral contributions from five different redshift bins, 
forz<0.5,0.5 <z< 1.5, 1.5 < z < 2.5, 2.5 < z < 4 andz > 4. 
This decomposition of the spectral contributions is also given 
for six different cosmogenic neutrino models. In all cases, 
Kotera et al. find that the neutrino fluxes are dominated by 
contributions from z > 1 .5, well into the asymptotic region of 
propagation time. Thus virtually all standard model cosmo- 
genic neutrinos have similar mean propagation times to Earth, 
(f v ) ~ 2 - 4 x 10 17 s, about a factor of 20-40 higher than that 
assumed by Cowsik et al. in their estimate. 



If a non-zero cosmogenic neutrino flux were observed, then 
solving equation [2] above, and requiring that the neutrino life- 
time should exceed the propagation time t v then gives an upper 
bound on a v : 



a v (£ v ) <4.0x 10 



-4 



5/3 -1/3 



^V.GeV 1 



(5) 



Inserting the typical value t v = 3 x 10 17 seconds, and convert- 
ing energy to EeV units (1 EeV = 10 9 GeV) which are more 
appropriate for the cosmogenic neutrino range 



oc v (£v) <6x 10 



-25 



£v V 5/3 
1 EeV J 



(6) 



which, as Cowsik et al note, is many orders of magnitude be- 
low the current best bound from IceCube observations of at- 
mospheric neutrinos. 

In fact this analysis is incomplete. Since the observed en- 
ergy at Earth E v /, s is redshifted relative to the energy at the 
source, in fact E v o b s — E v , src / (1 +z). Given that the domi- 
nant redshift range for cosmogenic neutrinos is 1 < z < 4, the 
bounds on oc v are factors of 3-15 lower; we can approximate 
this as 



a v (£v) <2x 10 



-25 



(z)E v 



obs 



1 EeV 



-5/3 



(7) 



where (z) is the weighted-mean z of the cosmogenic neutri- 
nos for a given source, and 1 < (z) < 4 is the valid range. This 
somewhat overestimates the bound since the propagation does 
not all occur at high-z; a more precise estimate requires inte- 
grating the energy dependence as a function of z, as we will 
show below, but equation[7]is a reasonable approximation. 
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FIG. 2: The attenuation coefficient for a narrow range of source neu- 
trino energies (eg. un-redshifted energies) for the given value of oc v , 
using the high-lambda cosmology and equation|9] 



Upper limits on fluxes of cosmogenic neutrinos do not of 
course uniquely point to new physics in order to suppress the 
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flux, unless the limits begin to strongly contradict other ob- 
servations. For example, because the GZK process relies on 
well-established particle physics with lab-frame energies in 
the GeV range, along with the extremely well-measured cos- 
mic microwave background radiation, the observed UHECR 
fluxes guarantee an associated neutrino flux, if the primary 
UHECRs are of a light composition, and the UHECR sources 
are distributed like other sources in the universe. A predom- 
inantly heavy (eg. iron) composition for the UHECRs, such 
as that suggested by Auger Observatory measurements 11281 . 
will significantly suppress the cosmogenic neutrino fluxes, 
but only if the the same UHECR composition measured for 
sources within the local universe applies to cosmological 
sources at z > 1.5. There is also currently some inconsistency 
in assuming a heavy composition for the primary UHECR, 
since the observed GZK cutoff would require significant tun- 
ing of the source energy spectra in order to match the cutoff, 
which would arise naturally for a light composition. 

Despite these uncertainties, the infinitesimal level of LIV in 
the neutrino sector required to effectively kill the cosmogenic 
flux suggests that we should take this possibility seriously 
and consider its implications. Within the last several years, 
a number of experiments, including the Radio Ice Cherenkov 
Experiment(RICE) E3, IceCube ED, and ANITA |[3Tll32l . 
have finally begun to constrain fluxes of cosmogenic neutri- 
nos. Although the total model space is not yet overly re- 
stricted by these constraints, the models that are ruled out, 
termed strong source evolutionary models, are still consistent 
with current UHECR data. These models do come into ten- 
sion with indirect bounds derived by analyzing extragalactic 
gamma-ray data from the Fermi satellite [33 , 41], but these de- 
pend on details of the secondary cascade process and UHECR 
source characteristics and these fluxes are thus not uniquely 
excluded fl2l . We therefore take these model bounds as mo- 
tivation to investigate what degree of LIV is required to sup- 
press these model fluxes below the current limits, and what 
other observable implications this might have. 

For A' neutrinos of energy E y src propagating from the cos- 
mogenic source to Earth, the LIV-induced attenuation with 
propagation time is given by dN/dt = —XN with the decay 
length X — l/x v . Then 



dN_dNdt_ r\(z)N 
dz 



dt dz T v (E v (z)) 
Substituting for x v , rearranging and integrating both sides 



N(z) = N exp 



«v 



dz'T](z') 



1 



(8) 



(9) 



Figure [2] shows an example of the attenuation coefficient that 
arises from this propagation equation, for neutrino source en- 
ergies over a quite small range, from 3-7 EeV, and a value 
of a v = 8 x 10~ 26 (chosen at random within the range of in- 
terest here) for the LIV parameter. Here the redshift is that 
seen from the source looking out toward Earth, using the high- 
lambda model. It is evident that in each case the attenuation 
has saturated by z ~ 1 relative to the source, which means 
that any typical cosmogenic neutrino model will undergo the 



maximal attenuation of its beam toward Earth. We empha- 
size the very strong energy dependence displayed here, with 
attenuation for a given value of a v increasing by two orders 
of magnitude over an octave of energy. This will tend to have 
a "brick-wall" impact on neutrino model fluxes, cutting them 
off very rapidly above some onset energy, as we show below. 




log io (E v /GeV) 



FIG. 3: Example of LIV effects on a strong-source-model cosmo- 
genic neutrino spectrum. Black line: unattenuated cosmogenic neu- 
trino spectrum generated by our modified CRpropa code, here with 
parameters chosen to give a spectral shape commensurate with that 
of typical strong-source cosmogenic neutrino models. The different 
lines indicated in the legend show the value of log 10 (a v that yields 
the modified spectrum shown. In each case both the brick-wall cut- 
off at a given energy, and an enhancement or pile-up effect just below 
that energy, appear in the LIV-modified spectra. 

We note also that the fact that there is a neutrino in the fi- 
nal state of this interaction, carrying on average 1/4 of the 
energy, means that the spectrum will not be simply attenu- 
ated above the cutoff energy, but that there should also be a 
pile-up of lower energy neutrinos just at the edge of the effec- 
tive LIV cutoff for a given values of a v . Both the brick-wall 
and pile-up effects are confirmed by an analysis from Bi et al. 
(2011) [36 1, where much lower energy neutrinos are consid- 
ered, over galactic propagation scales. 

To determine the impact of both the energy-dependent at- 
tenuation and the potential final-state neutrino pileup on an 
UHE neutrino spectrum at earth, we use the publicly avail- 
able cosmic-ray propagation code CRPropa 2.0 [37] to sim- 
ulate one-dimensional proton trajectories and then propagate 
the resulting neutrinos ourselves. We simulate protons from 
cosmic monoenergetic point sources out to z=3 [38j. Us- 
ing the energy and redshift of each neutrino produced along 
the trajectory, we allow it to interact via pair bremsstrahlung 
along its path to earth. We subtract the cumulative attenua- 
tion factors determined by equation[9]from unity to obtain the 
energy-dependent probability density functions used to find 
the redshift of each neutrino interaction. Multiple interactions 
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are allowed from each initial neutrino. Using the resulting 
collection of spectra from monoenergetic point sources, we 
follow the prescription for typical strong-source evolutionary 
models to calculate the neutrino spectra for an arbitrary in- 
jection spectrum and redshift evolution. We do not yet in- 
corporate possible energy-dependence of the LIV parameter 
a v , although such dependence may be expected from quan- 
tum gravity considerations [2]; in any case such effects would 
increase the sensitivity of our results to the level of LIV. 
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FIG. 4: ANITA-II and RICE limits and three strong-source-evolution 
model fluxes addressed here. 

For the final-state outgoing neutrino, we assume the pair 
bremsstrahlung process as a mean inelasticity such that the 
final state neutrino's energy is 1/4 of the initial energy. We 
would normally integrate over the inelasticity distribution, but 
such distributions are not given in the current models. In gen- 
eral we have found that, at least for neutrino deep-inelastic 
scattering inelasticity distributions, using the mean value as 
a proxy for a full integral is an acceptable approximation to 
first order. Fig. [3] shows an example of these results, here ap- 
plied to a cosmogenic neutrino model at varying levels of a v 
ranging from 10~ 27 to 10~ 21 . Both the brick-wall and pile-up 
effects are evident, and the turnover energy of the brick-wall 
in each case reflects the value of a v present. 

Returning to the strong-source-evolution models currently 
constrained by ANITA & RICE, Fig. [4] shows the fluxes for 
three of these models, all of which are constrained above the 
90% confidence level; also plotted are the ANITA and RICE 
limit curves. Energy flux (or intensity) units (particle energy 
x flux) are used here to avoid the spectral distortions that 
accompany less physically-motivated units such as E 2 F(E) 
which are also often used in reporting neutrino limits. In the 
normalization used here, models which just match the differ- 
ential limit curve for a decade of energy produce about 2 de- 
tectable events; in these three cases, the expected number was 



between 4-6 events for ANITA for all three models. For RICE, 
we estimate that the highest model of Barger et al. would pro- 
duce of order 12 events; the higher-energy models will pro- 
duce 2-3 events in RICE. 

TABLE I: First column: expected numbers of events N v for ANITA- 
II & RICE from three strong-source-evolution cosmogenic neutrino 
models; these models are all excluded at > 90% confidence from 
these experiments. Second column: values for the level of Lorentz 
invariance violation that would lead to the non-observation of these 
models. The high-lambda cosmological model was employed for this 
calculation. 



Model & references 


ANITA/RICE 


lower limit on 




predicted N v 


LrVctv 


Barger et al. 2006 [27] 






ANITA-II 


3.5 


2 x 10~ 27 


RICE2011 


12 


2 x 10~ 25 


Berezinsky 2005 (33) 






ANITA-II 


5.1 


3 x 10~ 28 


RICE2011 


3.4 


2 x 10~ 27 


Kalashev et al. 2002 (34l 






ANITA-II 


5.6 


2 x 10~ 28 


RICE2011 


2.9 


10 -28 



The event totals produced by an LIV-modified spectrum are 
determined for ANITA using ANITA event-sensitivity inte- 
grals based on our system Monte-carlos [39 1; estimating these 
event totals for RICE was outside the scope of our effort, but it 
is evident based on Fig.|4]that the constraints from RICE from 
the Barger et al. model could be a factor of three better for 
that case at least. Table[l]gives the hypothetical lower limits on 
0C V for each of the three cosmogenic neutrino models - lower 
limits here because they are the minimum values required to 
suppress these fluxes below the experimetal detection levels 
for RICE and ANITA-II. Not surprisingly the model with the 
highest energy predictions has the highest sensitivity to the 
LIV parameter. In all cases, even extremely small values of 
LIV, far below existing limits in the flavor-blind case, would 
be adequate to make these cosmogenic fluxes essentially un- 
detectable, accounting for the non-observation of these fluxes. 

It is an interesting coincidence that the level of LIV required 
to produce a flavor-blind superluminal effect that would effec- 
tively suppress cosmogenic fluxes is comparable to the current 
best limits on velocity differences between mass eigenstates, 
as noted in the introduction. Of course there are other poten- 
tial causes for suppression of the cosmogenic neutrino flux, 
and thus detection of the cosmogenic neutrinos is clearly a 
more interesting outcome for both neutrino astrophysics and 
for the resulting constraints on LIV phenomena. We can thus 
estimate the strength of the upper bound on LIV for a hypo- 
thetical detection of cosmogenic neutrinos for the sensitivity 
predicted for ANITA-III, due to fly in 2013, under the assump- 
tion of near-zero background, as has been the case for the 
first two ANITA flights. We have developed full Monte Carlo 
simulations of the sensitivity of the ANITA-III payload cur- 
rently in fabrication, and based on the planned improvements 



6 



we estimate about a factor of three increase in total neutrino 
event rate for a typical 30-day flight compared to ANITA-II. 
Under this outcome, with a clear detection of several neu- 
trino events, a conservative constraint on the LIV parameter 
from ANITA-HI would be a v ^ 10~ 26 , which is consistent 
with equation [7] above for (z) ~ 2 and £v,ofes ~ 3 EeV. Sim- 
ilar constraints would obtain for future detection of several 
unambiguous neutrino candidates in RICE, or the future ARA 
experiment now in early construction BUI , or for the recently 
fully-instrumented IceCube experiment. Thus a detection of 
cosmogenic neutrinos will move the flavor-blind constraints 
to values very close to those that currently obtain between the 
mass eigenstates. 

As noted above, a specific prediction of such LIV effects is 
the very steep dropoff with energy once a threshold is reached. 
Thus, large-scale experiments such as IceCube [30], Askaryan 
Radio Array (ARA) [40], and others that have extended low- 
energy sensitivity to cosmogenic neutrino models as com- 
pared to ANITA, may be able to detect such sharp cutoffs in 
neutrino energy spectra, signaling the onset of LIV at these 
extremely small levels. Thus our results modify somewhat 
the conclusion of Cowsik et al.: a clear detection of cosmo- 
genic neutrinos with no unexpected spectral features would 
provide exceedingly tight constraints on LIV in the neutrino 
sector. However, if enough neutrinos were measured to in- 
dicate the very sharp energy-spectral cutoff predicted here, it 
would rather unambiguously indicate the presence of neutrino 
LIV; no other proposed process that we know of can pro- 
duce such a sharp spectral cutoff in UHE neutrinos. These 
results emphasize the need for a broad-spectral approach to- 
ward searches for the cosmogenic neutrinos; this particular 
effect will turn on first at the highest energies, and might not 



be apparent in a cosmogenic neutrino detector whose spectral 
reach did not extend beyond 10 EeV or more. 

We conclude by observing that, although a relatively large 
value of oc v could completely attenuate the cosmogenic neu- 
trino fluxes below any detectable level for any experiment, the 
results of this would not escape some other possible detec- 
tion channels. Specifically, the pair bremsstrahlung process 
will lead to electromagnetic cascades in intergalactic space, 
redistributing the neutrino energy into an observable gamma- 
ray background. For example, in many models, the photo- 
hadronic cascades that UHECRs undergo in the GZK process 
lead to rough equipartition in energy between EeV cosmo- 
genic neutrinos, and GeV cosmogenic gamma-rays. In the 
presence of complete LIV-induced attenuation of the cosmo- 
genic neutrinos, the cosmogenic gamma-ray fluxes would be 
doubled to first order. Thus constraints on such a process 
can be obtained through careful analysis of the extragalac- 
tic gamma-ray background, such as that measured by Fermi, 
in direct analogy to what has already been done to constrain 
the integrated intensity of the cosmogenic neutrino produc- 
tion 11411 . Such constraints may in fact be more robust than 
those on GZK-interaction cascades, since LIV-induced neu- 
trino cascades are much less influenced by uncertainties in 
magnetic fields in the UHECR source environment. 
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